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Abstract 

We investigate the quantum mechanical osciUations of neutrinos prop- 
, agating in weak gravitational field. The correction to the result in the flat 

' space-time is derived. 



Q ■ Recently, the neutrino oscillation phase due to the presence of gravity is 

D , much discussed Q, |2|, Q. If the neutrinos are massive particles and mixed, 

Q I the neutrino oscillation between different flavor states occurs. Suppose that 

■ neutrinos are created as weak flavor eigenstate, say, at f^, and propagate to fe. 
. The state is a linear superposition of mass eigenstates and each phase of mass 
I eigenstate evolves in different way. As a result, the mixing phase angle for the 

J> ■ relativistic neutrinos propagating in a flat space is given by (fo = Arn^L/iihE), 

^ . where E is energy, L = \ fB — r^l, and Am^ = mf — m\. (See e.g., [^.) 

\ Gravitational effect was not seriously studied so far. Gravity can be elimi- 

CN ' nated by choosing appropriate inertial frame locally. The propagating distance 

■ of the neutrinos is so long in some case, that gravitational effect may become 
Q"^ \ important. Ahluwalia and Burgard |]|] considered the gravitational effect on the 

I neutrino oscillation. They showed that the external weak gravitational field of 

O^' a star with mass M adds a new contribution to the phase difference, denoted by 

. = —A.m?L{(j))/{AhE), where {(p) is defined by the average of gravitational 

potential over the semi-classical path, i.e., {4>) = —{J^^dLGM/r)/L. This 

' gravitationally induced phase can be estimated as (fc = ~{<P)fo- The phase 

^ . becomes to the extent of roughly 20% of ipo, near neutron stars. They suggested 

I that the new oscillation phase may be significant effect on the supernova explo- 

- - ' sions, since the extremely large fluxes of neutrinos are produced with different 

energies corresponding to the flavor states. 

The gravitationally induced oscillation phase may have the important astro- 
physical consequences. However, their derivation and even the definitions such 
as energy were not clear in their original paper. Bhattacharya, Habib and Mot- 
tola critically re-examined the quantum mechanical phase mixing. They 
calculated the phase difference for radially propagating particles, and found 
that the term of Avn? L{4)) / {UE) ~ GMAm^ log(rs/r^)/(7i_E) is canceled out. 
They showed that the possible gravitational effect appears at the higher order, 
Am^/E'^, and that the phase difference for radially propagating particles is 
GMAm^ log (r^/r^) /{4hE^). Numerically its magnitude is equal to ~ 10~^ 
(M/Mq) {Am^/eV^) (£'/MeV)~^ log(rB/r^), which is completely negligible 
in typical astrophysical applications. Therefore, the conclusion of Ahluwalia 
and Burgard [|^ seems to be incorrect for radially propagating case. Natural 
question is what happens in more general case. Does the term of GM Am? / E 
always disappear? 
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Furthermore, several authors have discussed the possibihty of the violation 
of equivalence principle in the neutrino oscillation, e.g.,[^], and reference 
therein. If the universality of the gravitational couplings to different flavors 
breaks down, additional phase difference appears. It is, therefore, an important 
matter to understand the modification of the neutrino oscillation phase due to 
the presence of gravity, even within the equivalence principle. This question 
must be settled first. In this paper, we will consider the quantum mechanical 
phase mixing of the neutrinos propagating in the weak gravitational field in 
detail. 

Let us consider a weak flavor eigenstate (a = e, ^u, or r), created at va 
with a certain energy. The value can be specified by the energy at infinity E. 
The flavor eigenstate is a linear superposition of mass eigenstates, which are 
represented by 

Wa) UaaWa), (1) 



where Uaa is unitary mixing matrix. The different mass eigenstates propagate 
with different velocity. The classical paths deviate reciprocally in the gravita- 
tional field. We use wave packet formalism (0, P) to examine the evolution 
of mass eigenstates. We assume that each mass eigenstate is created with the 
same energy as the Gaussian wave form in momentum space. The direction 
of the momentum is the same, but the mean value is different due to different 
mass. The wave packet for mass eigenstate, a can be written in coordinate 
space as 



iS, 



\^ ^ |2 



^aiS,t) = (2w^)-^/4exp (^^j exp ^-L_^j , (2) 

where S is phase function and should satisfy g^'^ S,fj_ Sy + m'^ = in geometrical 
optics limit. The size of the wave packet, a is assumed to be much smaller than 
the curvature radius of the external gravitational field. The center of the wave 
packet, Za{t) is determined by classical geodesic. 

The quantum mechanical transition probability from [va) to is 



aa I 



(2w2)-3/2 J2 UpaKaUphKb CXp [-{Sa - St) F, (3) 
ab ^'^ ^ 



where 



F{x,t) = exp - Zap + |f - 

= exp - ^(i*a + exp - . (4) 

The function F has a peak at {za + zi))/2, so that the phase factor can be 
evaluated by taking the dominant contribution in the stationary point. The 
factor exp(— — Zf,p/(8cj^)) represents that the coherence will be lost if two 
orbits separate much more than the size of the wave packet. The unitary matrix 



2 



for the mixing between two generations is parameterized by 6. If the deviation 
of the orbit is negligible, then we have 



(5) 



where (/? is the phase angle between two states, (p = 
In the opposite case, we have time-averaged one, 

P,_^ = ^sin2(20). 



-AS/i2n) = {S2-Si)/{2h). 



(6) 



We now estimate the split of the orbit, zi— Z2, and the phase difference AS" = 
—2hip, between two mass eigenstates. The classical trajectory of a particle with 
mass m around non-rotating star with mass M can be described by the energy 
at infinity E and angular momentum L. The orbital plane may be chosen as 
the equatorial plane of the spherical coordinate. Then the trajectory leaving at 
t = from (r, x) = {ta, 0) can be written as 



X 



rA 



E 



rA 



2GM\-^ 



dr, 



where 



B = E^ 



1 




(7) 
(8) 

(9) 



We will consider unbounded orbit, i.e., E ^ m. It is convenient to use the 
turning point, ro, of the radial direction to calculate the above integrals. Elim- 
inating L and expanding by G, we have 



X 



sm h 



GM 2r + ro 



1 



ro qGM 



rro 



r + ro 1 



rro 



rA 



(10) 




ro , {2-3q)GM 

+ --. log 

1 - a 



r + 



rA 
(11) 

where q = m? /E"^, and [/(r)]^!^ = /(r) — /{ta)- The first terms in eqs. (p^) 
and (11) represent straight lines in the spherical coordinate. The first order 
terms in G are well known effects on the propagation of massless particles. The 
second term in ([lO| ) is relevant to deflection of light, and the logarithmic term 
in ( pT[ ) is relevant to the time delay in the radar echo experiment. (See, e.g., 
§•) 

The phase S is given by the four momentum, conjugate to as. 



S 



-Edt+pidx" 



(12) 
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Both mass eigenstate are assumed to have the same conserved energy E, so 
that the first term / Edt is canceled out. Q We only calculate the contribution 
from the second term. The phase function can be evaluated for the classical 
path of the particle with mass m and energy E as 



X 



- + GM \ j — - log 




r + ro 



Expanding by m?/E'^ for the relativistic particles, we have 



5 ~ ^ 



r^ + GM 21og 



m 
'2E 

4 

m 



r2 - - GM I 2 log 




+ 



It is clear that the logarithmic term never appears in order GMmP'/E, but m 
order GMm^ / E'^ . However, there is still another term in order GMm? / E. We 
only consider the gravitational effect up to GMm? / E, and neglect the higher 
order terms hereafter. The particle with m? it Am2/2 arrives at (r it Ar, x i 
Ax/2) at the same coordinate time t. The split of the orbit can be calculated 
from eqs. ( p^ ) and (pT|). The phase difference A 5 can be calculated from 
eq.(|l^) at r. We explicitly show the differences in the path and phase due to 
mass square difference for the relativistic particles emitted to the radial and 
azimuthal directions, i.e., X = and r = initially. The turning points for 
these orbits can be chosen as tq = 0, and ro = r^, respectively. 



{1) radial orbits 

The orbital deviation and phase difference in the radial direction are given 

by, 

Am2 2GM] ^ 

Ar ~ — < 1 } L, (15) 

r J 



2^2 



A5~ 



Am^ 

where L is the path length in flat space, L = 



-L, 



(16) 



^ If the mixed state is a superposition of states with different energy, there is additional 
phase difference, AS = —AEt = —lS.m?t/{2E). We have the logarithmic term in AS by 
eq.(^f|), in converting from the coordinate time to the propagating distance. We will not 
consider such a term from the transformation between t and r hereafter. 



{2)transverse orbits 

For the particles emitted transversely, we have 



Am2 r 4GM 3GM1 L , , 

Ar=.--— 1 + — -L, 17 

2i/^ 1^ r + r^i r J r 

Am2 3GM GM 2GM 2GMr\ 



^^X^-ITT^ 1 + + — —\—^^ (18) 



where L = yr — rj^ in this case. The differences increase with the prop- 
agating distance L, as expected. Gravitational corrections are shown in the 
braces {}. The coherent wave packet with the initial size a will survive un- 
til L ~ E'^a I Arr? . Gravitationally induced oscillation phase appears in order 
GAm? / E except the purely radial motion, in which the gravitational correc- 
tion is accidentally canceled as seen in eq. dH). The magnitude of the oscillation 
phase can be written as (pc ~ {<P)fQ, from the dimensional argument. How- 
ever, the factor (0) is not the average of the gravitational potential over the 
semi-classical path. Rather, the correction can be regarded as the modification 
of the propagating distance. 
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